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The high density (HD) behaviour of nuclear symmetry energy E sym (p) is very impor- 
tant for understanding many interesting astrophysical phenomena, but it also subjects to 
the worst uncertainty among all properties of dense nuclear matter |T|]. To study the HD 
behaviour of E sym (p) has been a longstanding goal of extensive research with various mi- 
croscopic and/or phenomenological models over the last few decades, e.g., 0, for a recent 
review, see, e.g. The model predictions are extremely diverse and often contradictory. 
The results can be roughly classified into two groups, i.e., a group where the E sym (p) rises 
and one in which it falls with the increasing density. For instance, within the Hartree-Fock 
approach using about 25 Skyrme and Gony effective interactions that have been widely 
used successfully in studying saturation properties of symmetric nuclear matter and nuclear 
structures near the /3 stability valley, the calculated symmetry energies were found to fall ap- 
proximately equally into the two groups |||| . The density dependence of nuclear symmetry 
energy, especially at high densities, has many profound consequences for various studies in 
astrophysics [§,0. In particular, an increasing E sym (p) leads to a relatively more proton- rich 
neutron star whereas a decreasing one would make the neutron star a pure neutron matter 
at high densities. Consequently, the chemical composition and cooling mechanisms of pro- 
toneutron stars [|,[| , critical densities for Kaon condensations in dense stellar matter |TU| , |IT 



mass-radius correlations |12| , |13| as well as the possibility of a mixed quark-hadron phase |14 
in neutron stars will all be rather different. The fundamental cause of the extremely uncer- 
tain HD behaviour of E sym (p) is the complete lack of terrestrial laboratory data to constrain 
directly the model predictions. In this Letter, high energy heavy-ion collisions are proposed 
as a novel means to constrain stringently the HD behaviour of E sym (p). The upper bound 
of HD symmetry energy effects on high energy heavy-ion collisions is explored. It is shown 
for the first time that the neutron/proton ratio of HD nuclear matter formed in high energy 
heavy- ion collisions is uniquely determined by the HD behaviour of E sym (p). Experimental 
signatures of the HD behaviour of E sym (p) in two sensitive observables, i.e., the tt~ /h + ratio 
and the neutron-proton differential collective flow, are also investigated. 

We use here the parabolic approximation of the equation of state (EOS) for isospin 



asymmetric nuclear matter |3j 

e(p,5) = e(p,0) + E sym {p)5 2 , (1) 

where e(p, 0) is the energy per nucleon in symmetric nuclear matter, and 5 = (p n — p P )/(p n + 
p p ) is the isospin asymmetry in terms of the neutron (p n ) and proton (p p ) density, respec- 
tively. For the isoscalar part we use the simplest, momentum-independent parameterization 

e(p, 0) = a -u + j^u' + \e F u 2 '\ (2) 

where u = p/po is the reduced density and e° F = 36 MeV is the Fermi energy. The parameters 
a = —358.1 MeV, b = 304.8 MeV and a = 7/6 are determined by saturation properties and 
a compressibility = 201 MeV of isospin symmetric nuclear matter. Our conclusions in 
this work are independent of the particular form of e(p, 0). The symmetry energy is |l]|| 

E sym (p) = e(p, 1) - e(p, 0) = -E kin (p, 0) + V 2 (p), (3) 

where Eki n (p,0) is the kinetic energy per nucleon in symmetric nuclear matter and V^p) 
is the deviation of the interaction energy of pure neutron matter from that of symmetric 
nuclear matter. The E sym (p) becomes negative if the condition V 2 (p) < — ^Eki n (p, 0) is 
reached at high densities. Thus a pure neutron matter could become most stable, leading 
to the isospin separation instability in HD neutron-rich matter. Consequently, pure neutron 
domains or neutron bubbles surrounding isolated protons can be formed in neutron stars 
|J. To represent the two groups of model predictions, we use for the E sym (p) 

E sym(p) = E sym (p )u and E sym (p) = E sym (p )u ■ c _ (4) 

where E sym (p ) = 30 MeV is the symmetry energy at normal nuclear matter density, and 
u c = Pel Po is the reduced critical density at which the E b sym (p) crosses zero and becomes 
negative at higher densities. The predicted value of u c ranges from about 2.7 (Hartree- 
Fock with the Skyrme interaction Sp 0) to 9 (variational many-body approach with the 
UV14+UVII interaction [§].)• These two forms of the E sym (p) (insert in the lower window) 
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together with the corresponding EOS for isospin asymmetric nuclear matter are shown in 
Fig. 1. With the linearly increasing E^ ym (p), the EOS (upper window) becomes stiffer with 
the increasing 5. The isospin symmetric nuclear matter remains to be the ground state at 
all densities. This is in stark contrast to the situation using the E b sym (p). The EOS obtained 
with the E b syrn (p) and u c = 3 (lower window), is softened instead of being stiffened by the 
increasing isospin asymmetry 8 at densities higher than 3po- At these high densities, the 
pure neutron matter and the isospin symmetric nuclear matter thus becomes the most stable 
and unstable form of HD nuclear matter, respectively. 

High energy heavy-ion collisions provide the only terrestrial situation where the HD 
neutron-rich matter can be formed. Moreover, fast radioactive heavy-ion beams to be avail- 
able at the planned Rare Isotope Accelerator (RIA) in the United States can make the HD 
nuclear matter even more neutron-rich. The HD behaviour of E sym (p) affects properties of 
the HD nuclear matter formed in high energy heavy-ion collisions. Moreover, it leads to 
interesting precursory phenomena in the experimental observables of high energy heavy-ion 
collisions. We investigate these effects and phenomena within an isospin-dependent hadronic 
transport model |TR[T(1. Evolutions of the phase-space distribution functions of nucleons, 
Delta resonances and pions with their explicit isospin degrees of freedom are solved nu- 
merically by using of the test-particle approach [pLTl, 18|1 - Both the isoscalar and isovector 
mean-filed potentials are derived consistently from the nuclear EOS given above. Isospin- 
dependent total and differential cross sections among all particles are taken either from 
the elementary particle scattering data or obtained by using the detailed balance. Explicitly 
isospin-dependent Pauli blockings for fermions are also employed. For a review of the model, 



we refer the reader to ref. 19 



Shown in Fig. 2 is the correlation between the isospin asymmetry 5 and the baryon 
density p/po at the instant of 20 fm/c in the reaction of 132 Sn + 124 Sn at a beam energy of 
400 MeV/nucleon and an impact parameter of 1 fm. This very reaction will be available at 
the RIA facility. The overall rise of S at very low densities is due to the neutron skins of the 
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colliding nuclei. One also notices that the E b sym (p) leads to a slightly higher compression due 
to its relatively softening effect on the nuclear EOS. Dramatic effects due to the different 
symmetry energies are clearly revealed especially at high densities. For a comparison, the 
corresponding correlation in neutron stars at (3 equilibrium is shown in the insert. As a 
reference, an estimate of the direct URCA limit of 8p (below which the fast cooling process 
can happen in protoneutron stars and thus results in also higher neutrino flux) is indicated 
with the dotted line f|. To a good approximation perfectly suitable for this study, a neutron 
star can be considered as consisting of nuetrons, protons and electrons. At (3 equilibrium 
the proton fraction xp = p p / (p n + p p ) in neutron stars is then determined by || 

ftc(37rV/3) 1/3 = 4£ sym (p)(l - 2x p ). (5) 

The equilibrium isospin asymmetry 5/3 = 1 — 2xp is therefore entirely determined by the 
E S ym(p)- With the E^ ym (p), the 5p is 1 for p/p > 3, indicating that the neutron star has 
become a pure neutron matter at these high densities. On the contrary, with the E^ ym (p), the 
neutron star becomes so proton-rich that the direct URCA process can happen at densities 
higher than about 2.3po- An astonishing similarity is seen in the resultant 5 — p correlations 
for the neutron star and the heavy-ion collision. In both cases, the symmetry energy E h sym {p) 
makes the HD nuclear matter more neutron-rich than the E® ym (p) and the effect grows with 
the increasing density. Of course, this is no surprise since the same nuclear EOS is at work 
in both cases. As shown in Fig. 1, in stark contrast to the situation with the E® ym (p), the 
decreasing E h sym (p) above po makes it energetically more favorable to have the denser region 
more neutron-rich. 

To further investigate the HD behaviour of E sym (p), the average n/p ratio of the HD 
region with p/po > 1 is examined as a function of time and beam energy in the left panels 
of Fig. 3. The effect on (^/p) p / po >i due to the different E sym (p) is seen to grow with both 
the reaction time and the beam energy. This is because of the higher densities reached with 
the more energetic beams and after longer times of compression. 

How to probe experimentally the HD behaviour of E sym (p) in high energy heavy- ion 
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collisions? Among the observables we have explored, the tt~/tt + ratio and the neutron- 
proton differential collective flow are found to be most promising. The former measures 
sensitively the neutron to proton ratio n/p of the HD nuclear matter and thus indirectly 
the HD behaviour of E sym (p). At beam energies below about 2 GeV/nucleon, pions are 
mostly produced through the decay of A(1232) resonances. The primordial tt^/tt + ratio is 
approximately quadratic in n/p according to the branching ratios of single pion production 
via A resonances in nucleon-nucleon collisions 

vr-A + = ^±1* « {n/pf. (6) 
hp 2 + np 

Pion reabsorptions and rescatterings are expected to complicate the above relationship. 
Nevertheless, very high sensitivity to the n/p ratio is retained in the final tt~/tt + ratio as 
indicated in the experimental data of high energy heavy- ion collisions |]20[| . Shown in the 
right panels of Fig. 3 are the (n~ /n + )uk e ratio 

(* A )«* = 7r+ + A++ + iA+ (7) 

as a function of time. This ratio naturally becomes the final 7r~ /tt + ratio when the reaction 
time t is much longer than the lifetime of the delta resonance ta- The (7T~ /7r + )n^ e ratio is 
rather high in the early stage of the reaction because of the large numbers of neutron-neutron 
scatterings near the surfaces of the colliding nuclei. By comparing the two sides in Fig. 3, 
it is seen that a variation of about 30% in the {n/p) p / Po >i ratio due to the different E sym (p) 
results in about 15% change in the final 7t~/tt + ratio. The later thus has an appreciable 
response factor of about 0.5 to the variation of HD n/p ratio, and it is approximately beam 
energy independent. 

The neutron-proton differential collective flow is measured by |^TJ 

Fnp(y) = Vxji, (8) 

where N(y) is the total number of free nucleons at the rapidity y, p Xi is the transverse 
momentum of particle % in the reaction plane, and Tj is +1 and —1 for neutrons and protons, 

6 



respectively. The F np (y) combines constructively the in-plane transverse momenta generated 
by the isovector potentials while reducing significantly influences of the isoscalar potentials 
of both neutrons and protons. Thus, it can reveal more directly the HD behaviour of E sym (p) 
in high energy heavy-ion collisions. A typical result for the 132 Sn+ 124 Sn reaction is shown in 
the upper window of Fig. 4. A clear signature of the HD behaviour of E sym (p) appears at both 
forward and backward rapidities. To characterize the effect, the slope dF np / d{y cm / ybeam) at 
mid-rapidity is shown as a function of beam energy in the lower window. A striking difference 
of about a factor of 2 exists in the reactions at Eb eam > 200 MeV/nucleon. This large effect 
can be very easily observed by using available detectors at several heavy-ion facilities in the 
world |2^j2^]. Compared with the 7t~/tt + ratio, the neutron-proton differential collective 
flow is more directly affected by and is thus also a more sensitive probe of the HD behaviour 
of E sym (p). 

In conclusion, the HD behaviour of E sym (p) has been puzzling physicists for decades. 
In this work, high energy heavy-ion collisions are proposed as a novel means to solve this 
longstanding problem. The upper bound of HD symmetry energy effects on high energy 
heavy-ion collisions is explored. For the first time, it is shown that the isospin asymmetry 
of HD nuclear matter formed in high energy heavy-ion collisions is uniquely determined by 
the HD behaviour of E sym (p). Both the n~ /n + ratio and the neutron-proton differential 
collective flow can serve as sensitive probes. Measurements of these observables will provide 
the first terrestrial data to constrain stringently the HD behaviour of nuclear symmetry 
energy. 

This work was supported in part by the National Science Foundation Grant No. PHY- 
0088934 and Arkansas Science and Technology Authority Grant No. 00-B-14. 
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FIGURES 




FIG. 1. The upper (lower) window is the equation of state of isospin-asymmetric nuclear matter 
corresponding to the nuclear symmetry energy parameterization Eg ym (E b sym ) shown in the insert 
of the lower window. 
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FIG. 2. The isospin asymmetry-density correlation at t=20 fm/c over the whole space in the 
132 Sn + 124 Sn reaction with the nuclear symmetry energy E a sym and E b sym , respectively. The 
corresponding correlation in neutron stars is shown in the insert. 
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FIG. 3. Left panels: the neutron/proton ratio of nuclear matter with density higher than 
normal nuclear matter density as a function of time with the nuclear symmetry energy E^ ym and 
Eg , respectively. Right panels: (tt~ /^ + )uke ratios in the same reactions. 
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FIG. 4. Upper window: the neutron-proton differential collective flow in the Sn + Sn 
reaction with the nuclear symmetry energy E% ym and E b sym , respectively. Lower window: excitation 
function of the slope parameter of the differential flow for the 132 Sn + 124 Sn reaction. The solid 
(with the Eg ym ) and dashed (with the E b sym ) lines are drawn to guide the eye. 



13 



